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Introduction: Atmospheric sulfur species resulting 
from volcanic emissions impact the composition and 
chemistry of the atmosphere, impact the climate, and 
hence, the habitability of Mars and impact the mineral- 
ogy and composition of the surface of Mars. The geo- 
chemical/photochemical cycling of sulfur species be- 
tween the interior (via volcanism), the atmosphere (at- 
mospheric photochemical and chemical processes) and 
the deposition of sulfuric acid on the surface of Mars is 
an important, but as yet poorly understood geochemi- 
cal/photochemical cycle on Mars. There is no observa- 
tional evidence to indicate that Mars is volcanically 
active at the present time, however, there is strong evi- 
dence that volcanism was an important and widespread 
process on early Mars. The chemistry and photochemi- 
stry of sulfur species in the early and present atmos- 
phere of Mars will be assessed using a one-dimensional 
photochemical model. Since it is generally assumed that 
the atmosphere of early Mars was significantly denser 
than the present 6-millibar atmosphere, photochemical 
calculations were performed for the present atmosphere 
and for the atmosphere of early Mars with assumed 
surface pressures of 60 and 350-millibars, where higher 
surface pressure resulted from enhanced atmospheric 
concentrations of carbon dioxide (C0 2 ). The following 
sections include the results of earlier modeling studies, 
a summary of the one-dimensional photochemical mod- 
el used in this study, a summary of the photochemistry 
and chemistry of sulfur species in the atmosphere of 
Mars and some of the results of the calculations. 

Earlier Studies: Modeling of early Mars suggests 
a warm, wet climate to explain the fluvial geomorphol- 
ogy seen on the surface of Mars [1]. To increase the 
atmospheric surface temperature, models of early Mars 
suggest high levels of atmospheric C0 2 [1,2]. Subse- 
quent analysis of C0 2 cloud formation causing surface 
cooling has caused inconsistencies in the C0 2 warming 
wet model’s approach, as well as the observed lack of 
carbonates found on the Martian surface [3,4,5]. 

Although differing in the chemistry, several scien- 
tists, beginning in 1997 with Yung et al. [6] proposed 
atmospheric sulfur dioxide (S0 2 ) as a potential solution 
for the early warm and wet climate of Mars. Yung et al. 
(6) proposed that with sufficient S0 2 present in the 
middle atmosphere,, C0 2 clouds would be abated and 
more surface warming could take place with higher 
C0 2 concentrations. 

Subsequent theories on the introduction of S0 2 ad- 
dress the additional question of sequestered carbonates. 
Acidic water from precipitation of H 2 S0 4 due to vol- 


canic release of S0 2 and H 2 S could halt carbonate for- 
mation and in turn lead to sulfur species formation as 
found presently in Martian surface data [5, 7]. Bullock 
and Moore [4] suggest acidic water from this precipita- 
tion would limit carbonate formation and lead to sulfate 
species precipitation as seen for example at West Can- 
dor Chasma [8], Meridiani Planum, and various regions 
in Valles Marineris [9] on Mars. SO? would help sus- 
tain a thick CO? atmosphere as described by Yung et al. 
[6] as limited carbonate formation would not sequester 
the C0 2 out of the atmosphere. Bullock and Moore [4] 
suggests this change in the atmosphere but did not 
model the effects of SO? in such an atmosphere for via- 
bility of this constituent’s lifetimes in model atmos- 
pheres. Thus, here a one dimensional steady state model 
will be used to address the introduction of volcanic sul- 
fur into the early Martian atmosphere. 


The Photochemical Model of the Mars Atmos- 
phere: The photochemical model that will be used in 
this research will be the model reported by Nair et al. 
[ 1 0] updated with additional methane and sulfur chemi- 
stry from Summers et al. [11] and Wong et al. [12]. 
This photochemical model, originated at Caltech/JPL, 
solves the one-dimensional continuity equation for 
steady-state conditions for the species listed in Nair et 
al. (10). The equation. 
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is solved where n, are the concentrations, <[>, are the 
vertical diffusive fluxes, and P, and L; are the chemical 
loss and production terms for each indexed species. 


Since the steady state term leads to — = 0 , the conti- 
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nuity equation can be solved giving the indexed species 
vertical flux by the equation: 
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where D , is the molecular diffusion coefficient, K is the 
eddy diffusion coefficient, T is the temperature, a is 
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thermal diffusion factor (which is taken to be zero), and 
H is the scale height of the background atmosphere 
[ 10 ]. 

The 27 species used in Nair et al. (10) have been 
augmented by Summers et al. [11] with an additional 37 
species. The addition of these species required addi- 
tional reactions to be added to the model. Where possi- 
ble JPL photochemical data was used. For a complete 
discussion of parameter value settings used in the mod- 
el see Nair et al. [10] and references within. 

The water profile shows both Nair et al. (10) high 
and low water mixing ratios and the saturation water 
level for the current temperature range. Water is initial- 
ly held constant within the model and below saturation 
under 40 km as shown in Figure 1. Although reactions 
were updated, initial model runs show agreement with 
Nair et al. [10]. Thus, changes can be made to represent 
early Martian atmosphere from this working current 
Martian atmospheric model. 



Figure 1. Calculated trace species in the atmosphere of 
Mars (C0 2 = 6 mbar). 


The Photochemistry and Chemistry of Sulfur 
Species in the Atmosphere of Mars: Once comparable 
results were demonstrated sulfur boundary conditions 
could be added to the model for both S0 2 and FESO4. 
S0 2 is proposed to enter the early atmosphere through 
volcanic emission while Fl 2 S04 condenses out of the 
atmosphere once the sulfur dioxide has cycled through 
chemical reactions in a proposed sulfur cycle shown in 
Figure 2 [5]. Rates for influx and deposition of chemi- 
cals were estimated through two distinct methods. S0 2 
influx, as a constituent of volcanic emission of similar 
Earth volcanoes, is dependent upon the amount of vol- 
canic activity and the constituent make up of the ex- 
pelled gases. Estimates have been made of the amount 
of S0 2 that entered the early atmosphere of Mars 
through analysis of surface geology and chemical make 
up of gases expelled from volcanoes on the Earth and 
compared to Martian volcanoes and chemical analysis 
of surface rock. Levine and Summers (5) estimate S0 2 


emissions for Mars to be between 4-8 x 10 11 moles 
S/year. This implies that an averaged upward flux rates 
would be approximately 1 x 10 10 cm's' 1 during the ear- 
ly volcanic period on Mars. This boundary condition 
was added to the model. 



Figure 2. The photochemistry and chemistry of sulfur 
species in the atmosphere of Mars. A surface boundary 
condition provides influx of S0 2 into the atmosphere is 
estimated from early volcanic emissions on Mars. 

With the introduction of S 0 2 into the atmosphere it 
is crucial to set a condensation rate for H 2 S04. The 
H 2 S04 must have an exit path in the atmosphere as well 
as a deposition rate on the surface of Mars. A pseudo 
reaction rate is found for FESO4 condensation and the 
surface boundary condition of FESO4 is fixed to be 
lxltr'cm/s, in line with this downward transport and 
precipitation out of the atmosphere and onto the sur- 
face. Models were then rerun to see affects of the intro- 
duction of sulfur species on the model. 

With the introduction of sulfur into the model the 
reaction to changes in atmospheric water should be 
examined. Three water scenarios will be used. Starting 
with the water profile used in the present day Martian 
model presented by Nair et al. [10] and augmented with 
water saturation curves, a new water profile is generat- 
ed for each of the three temperatures curves from Kast- 
ing [3] for each partial pressure of C0 2 corresponding 
to 6, 60 and 350-millibar atmospheres. 

All water vapor added to the model was limited by 
the water saturation curves for the given temperatures 
present at each elevation. This constraint is initially 
necessary so that no water clouds would be formed and 
therefore requiring additional changes within the at- 
mospheric model. The water profile above approx- 
imately 20 km is kept constant as the minimum water 
saturation seen for each model for simplicity. Sulfur 
species are not highly sensitive to the changes made to 
the water profiles with saturations limits. Now that wa- 
ter has been seen not to influence sulfur species we 
need to examine how the differing C0 2 levels affect the 
levels of sulfur species in the models. For these scena- 
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rios the water profiles will be reset to the lx Nair et al. 
(11) model augmented with the water saturation curve. 
Initial results demonstrate that the surface lifetime of 
S0 2 increases with increasing C0 2 . The calculated SO? 
lifetime increases by three orders of magnitude from 
approximately half a day to a bit over 1.5 years while 
C0 2 increases by less than 2 orders of magnitude (from 
6 to 350-millibars of C0 2 ) showing sensitivity to the 
major atmospheric constituent of present day Mars. 

The surface influx value estimated for sulfur dio- 
xide is dependent upon volcanic activity that is not 
tightly bounded both because of the variability of vol- 
canic activity and the basis of Earth’s volcanic output 
for the SO? estimates, therefore additional model runs 
need to examine the sensitivity of the results to changes 
in the sulfur species themselves. Since the surface in- 
flux is estimated from volcanic emissions of sulfur dio- 
xide, S0 2 influx values will be changed to mimic the 
variability of volcanic activity during active periods. 
Holding the water profile as outlined above for each 
model, changes were made in the S0 2 upward flux 
boundary condition to examine how the sulfur cycle 
would be affected by changes to the upward flux of 
S0 2 . The upward flux values were changed by an order 
of magnitude larger and smaller than the estimated flux 
value used to this point in the model. Changes by an 
order magnitude in the surface flux of S0 2 leads to sim- 
ilar flux. Therefore, the S0 2 flux profile appears unaf- 
fected by the change in the amount of S0 2 emission. 
Finally, inspection of the H 2 S0 4 flux profiles also show 
similar flux profiles except for changes between 10-20 
km. 

With S0 2 lifetime in the atmosphere increasing 
from approximately one day to one year with enhanced 
concentrations of carbon dioxide, the sulfur dioxide 
reacts with oxygen species or breaks up due to photoly- 
sis creating new sulfur species such as SO and S0 3 . 
These species are short-lived and relatively insensitive 
to the increase in C0 2 . H 2 S0 4 will be produced after 
S0 2 react with photons and oxygen species. It then ac- 
cumulates in the atmosphere and precipitates out onto 
the planetary surface. 

Since H 2 S0 4 accumulates in the atmosphere, a 
pseudo-reaction rate is introduced into the model to 
remove the sulfuric acid out of the atmosphere and 
model the precipitation on to the surface. This reaction 
is entered into the model reactions as: H 2 S0 4 + M — » 
PROD where M is a non-specific molecule and PROD 
is an end product removed from further atmospheric 
analysis. Thus, H 2 S0 4 molecules are effectively remov- 
ing from the atmosphere as desired. In setting a value 
for the reaction rate one must mimic the rate of conden- 
sation of H 2 S0 4 out of the atmosphere. To estimate this 
value we start by estimating the condensation time con- 
stant for H 2 S0 4 , t. This condensation time constant is 
translated into a reaction rate for entry into the model 


by solving &= kn n. . for the k where x = -7- . The 
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initial value used as the initial reaction rate for H 2 S0 4 is 
5 x 10" 23 cm 3 s’ 1 . Condensation rate sensitivities were 
also run for this pseudo-reaction rate, changing the rate 
by an order of 10, larger and smaller. The range of pa- 
rameters used in the photochemical calculations in this 
paper is summarized in table 1 . 

Results of Photochemical Calculations: 

1. Enhanced C0 2 concentrations significantly enhance 
the lifetime of atmospheric sulfur dioxide. 

2. Water mixing ratio limited at saturation does not 
influence the introduction of sulfur species with wa- 
ter, likely volcanic emissions. More work needs to 
be done with water saturation and cloud formation 
for definitive result. 

3. Surface flux boundary conditions do not show great 
sensitivity to order of magnitude changes to the up- 
ward flux of S0 2 . This is important as upward flux 
could change by orders of magnitude since volcanic 
activity was not constant over geological time 
scales. 

4. H 2 S0 4 formed in the atmosphere (the predominant 
product of S0 2 ) condensed onto the surface with the 
introduction of S0 2 into a predominately C0 2 at- 
mosphere — this could address the apparent lack of 
carbonates on the surface of Mars and explain the 
sulfur minerals found on the surface. 

Conclusions: Atmospheric sulfur species resulting 
from volcanic emissions impact the composition and 
chemistiy of the atmosphere, impact the climate, and 
hence, impact the habitability of Mars and impact the 
mineralogy and composition of the surface of Mars. 
The geochemical and photochemical cycling of sulfur 
species between the interior (via volcanism), the atmos- 
phere (atmospheric photochemical and chemical 
processes) and the deposition of sulfuric acid on the 
surface of Mars is an important, but as yet poorly un- 
derstood cycle on Mars. Many questions about sulfur 
species and their relationship to the geology, geophys- 
ics, atmospheric composition, chemistry and climate 
and habitability of Mars may have to await the future 
exploration of Mars by human explorers [12, 13, 14, 
15, 16], 
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Table 1: Range of Atmospheric Parameters for Photochemical Calculations in This Paper. 


Model Parameters 

Description 

Reference 

co 2 

a. SP = 6mbar 

Similar to current surface pressure P s = 0.006 bar 

Kasting [3] 

b. SP = 60mbar 

lOx La. 

Kasting [3] 

c. SP = 350mbar 

Maximum pressure before C0 2 cloud formation in 
stratosphere P s = 0.35 bar 

Kasting [3] 

Temperature 

a. Nair high temp, curve 

Model present day temp. 

Nair et al [10] 

b. Kasting La. temp, curve 

Model curve augmented with appropriately fitted 
Nair curve above 60km 

Kasting [3] 

c. Kasting Lb. temp, curve 

Model curve augmented with appropriately fitted 
Nair curve above 60km 

Kasting [3] 

d. Kasting l.c. temp, curve 

Model curve showing warmest temperature with 
only C0 2 in atmosphere and without cloud 
formation augmented with appropriately fitted Nair 
curve above 40 km 

Kasting [3] 

h 2 o 

a. Nair/saturation curve 

Nair water curve limited in values by H 2 0 
saturation curve for temperature curves above and 
remaining constant above 20 km 

Nair et al [10] 

Boundary Conditions S0 2 

a. surface flux BC = 1.06x10 10 cm'V 1 

Estimated flux value from comparison to Earth 
shield volcano output 

Levine & 
Summers [5] 

b. surface flux BC =1.06xl0 u cm'V 1 

1 Ox estimated amount 

Levine & 
Summers [5] 

c. surface flux BC = 1.06xl0° 9 cm'V 1 

1/1 0 th estimated flux 

Levine & 
Summers [5] 

H 2 S0 4 Condensation Rate 

a. pseudo reaction rate to model conden- 
sation of H 2 S0 4 = 5x1 0" 20 cm" 3 s" 1 

Estimated using theory from Summers and Strobel 

Summers & 
Strobel [17] 

b. pseudo reaction rate = 5xl0" 22 cm" 3 s" 1 

lOx estimated value 

Summers & 
Strobel [17] 

c. pseudo reaction rate = 5xl0' 19 cm" 3 s" 1 

l/10x estimated value 

Summers & 
Strobel [17] 
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